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ABSTRACT 

We observed diffuse cosmic and atmospheric gamma rays at balloon altitudes with the Sub-MeV 
gamma-ray Imaging Loaded-on-balloon Experiment I (SMILE-I) as the first step toward a future all- 
sky survey with a high sensitivity. SMILE-I employed an electron-tracking Compton camera comprised 
of a gaseous electron tracker as a Conipton-scattering target and a scintillation camera as an absorber. 
The balloon carrying the SMILE-I detector was launched from the Sanriku Balloon Center of the 
Institute of Space and Astronautical Science/ Japan Aerospace Exploration Agency (ISAS/JAXA) on 
September 1, 2006, and the flight lasted for 6.8 hr, including level flight for 4.1 hr at an altitude of 
32-35 km. During the level flight, we successfully detected 420 downward gamma rays between 100 
keV and 1 MeV at zenith angles below 60 degrees. To obtain the flux of diffuse cosmic gamma rays, 
we first simulated their scattering in the atmosphere using Geant4, and for gamma rays detected at an 
atmospheric depth of 7.0 g cm^^, wc found that 50% and 21% of the gamma rays at energies of 150 keV 
and 1 MeV, respectively, were scattered in the atmosphere prior to reaching the detector. Moreover, 
by using Geant4 simulations and the QinctiQ atmospheric radiation model, we estimated that the 
detected events consisted of diffuse cosmic and atmospheric gamma rays (79%), secondary photons 
produced in the instrument through the interaction between cosmic rays and materials surrounding 
the detector (19%), and other particles (2%). The obtained growth curve was comparable to Ling's 
model, and the fluxes of diffuse cosmic and atmospheric gamma rays were consistent with the results 
of previous experiments. The expected detection sensitivity of a future SMILE experiment measuring 
gamma rays between 150 keV and 20 MeV was estimated from our SMILE-I results and was found to 
be ten times better than that of other experiments at around 1 MeV. 

Subject headings: balloons — instrumentation: detectors — diffuse radiation — gamma rays: obser- 
vations 
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1. INTRODUCTION 

Observations in the low-energy gamma-ray band 
from hundreds of keV to tens of MeV provide a unique 
window on such phenomena as gamma-ray lines from 
nuclear de-excitation produced by nucleos ynthesis 
supernovae , e.g . , ^^Ni, ^^C o, and 
1988; Chevali eij 119921: Ivon B allmoos 
20011: iCheng fc Romerol [2004i; iBoggsi i2006|); gamma- 
ray lines from long-lived isotopes spr ead throughout 
our ga laxy, such as ^^ Al and ^ °Fe (lOberlack et al.l 
[l99l IKnodlseder eJall [19991: IDiehl et al.l I2003n : 
electron -positron ann i hilatio n lines from the Galactic 
center (liPurcell et al.' '19971 IKnodlseder et al.l 120031: 
IWeidenspointner eF al. 20 0fl); neutron capture lines 
from solar flares 111 are fc Murphvl 12001 : and gamma- 
ray resonant absorption along lines of sight toward 
gamma-ray bright quas ars due to delta-r esonance or 
giant dipole resonance (jlvudin et al.l I2005D . Further- 
more, radiation from black hole accretion disks has 



been detected (iMcConnell et al1l2002| ). and radiation of 
neutral pions produced by ions accelerated in the strong 
gravitational potential of black holes i s also expected to 
be d etected (iMcConneU et al.l Il996t iMahadevan et al.l 
Il997t iBhattacharvva et al.ll20'03[ r We can also observe 
nonthermal processes in low-energy gamma-rays: syn- 
chrotron radiation or inverse Com pton-scattered gamma 
rays in gamma-ray p uls ars (jAharonian fc Atovan 
19981: [Thompson et al.l 119991: iSchonfelder et al 



200a iKuiper et al.l 1200 ID. 



(AGN; 



active galactic nuclei 
jR)ssati et al. 199? 



. .Urrv fc Padovanil Il99l ._ 

Ghisellini et al.lll998HKubo et al.|[T99 8: Chiabcrgc ct al 
20011) ■ and gamma ray bursts (GRB; jBr iggs et al . 199^ 
Paciesas et all 119991 : iKaneko et all |2006|)! In addition 
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there are diffuse galactic gamma rays probably produced 
by bremss trahlung and inve r se Compton scatter ing of 
electrons (iB oggs et al] 120001: [Strong et all I2000D . and 
diffuse extragalactic gamma rays, which are thought to 
be a combinati on of emissions from A GNs and Type 
la supernovae ([The et al.l 119931 : iWatanabe et al.l Il999| ) 
or the combined emission from the Comptonization 
process including nonth ermal electrons in accretion 
disk coronae of AGNs ([Inoue et al.l I2008D . In addi- 
tion, the measurement of photon polarization through 
low-energy gamma-ray observations from nonthermal 
processes, such as synchrotron radiation and Compton 
scattering, would present a pow er ful diagnostic too l 
for astrophysics (|Lei et al.l 119971: iForot et al.l I2008D . 
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However, observation of this low-energy gamma-ray 
band is difficult because of the following reasons: the 
dominant process in a detector is Compton scattering 
and large backgrounds of photons are produced in the 
hadronic process between cosmic rays and a satellite 
body. Therefore, MeV gamma-ray astronomy has not 
advanced in comparison with X-ray or other gamma-ray 
bands. In fact, COMPTEL onboard the Compton 
Gamma Ray Observatory (CGRO) discovered only ^ 30 
steady gamma-ray sources in the 0.75-30 MeV band 
(jSchonfcldcr ct al. 2003), whereas EGRET detected 
^ 270 sources pfartman et al.l[T999| ). and Fermi found 
1451 so urces during the first 11 months of the all-sky 
survey (|Abdo et all I2010D in the sub-GeV/GeV region 
above 100 MeV. 

COMPTEL, which was the first Compton telescope on- 
board a satellite, localized the direction of an incident 
gamma ray on an event circle superposed on the sky by 
measuring the direction of a Compton-scattered gamma 
ray and the energies of both the scattered gamma ray 
and a Compton-recoil electron. After the launch, the 
sensitivity of COMPTEL was found to be mainly deter- 
mined not by diffuse cosmic gamma rays but by 1) locally 
produced background gamma rays; 2) non-Compton- 
scattering events such as multiple-photon events due 
to hadronic interaction, bremsstrahlung, or annihila- 
tion of positrons; and 3) events due to other particles 
such as neutrons, even after effective background re- 
jection by time-of-fiight measurement between a Comp- 
ton scattering material and an absorber of t he scattered 
gamma rays (jWeidenspointner et al.l [200l . COMP- 
TEL's limited sensitivity made it desirable to create a 
new instrument with a higher sensitivity in the sub- 
MeV/MeV region. Therefore, we developed an electron- 
tracking Compton camera (ETCC) consisting of a three- 
dimensional tracker of Compton-recoil electrons and an 
absorber of Compton-scattered gamma rays by a new 
detection niethod with powerful back ground rejection 
(jOrito et al.ll2003l:lTanimori et al.l200l . Compared with 
COMPTEL, the ETCC can restrict the direction of an 
incident gamma ray to a reduced arc on the Compton 
circle by measuring the three-dimensional track of the 
Compton-recoil electron in addition to the energy of the 
Compton-recoil electron and the energy and direction of 
the Compton-scattered gamma ray. Moreover, the angle 
between the direction of the scattered gamma ray and 
that of the Compton-recoil electron can be used for pow- 
erful background rejection by checking the consistency 
between the measured and kinematically calculated val- 
ues. The ETCC can also detect the polarization of an 
incident gamma ray, because the azimuth angle distribu- 
tion of Compton scattering has an asymmetry for polar- 
ized gamma rays. 

Because we had already confirmed the detection 
rinciple of t he ET CC by groun d-based experiments 
Takada et all 120051: iHattori et all l2007t) . we initiated 
an experiment to observe sub-MeV gamma rays from 
celestial objects by using a balloon-borne camera, the 
Sub-MeV gamma-ray Imaging Loaded-on-balloon Exper- 
iment (SMILE). For the first flight of SMILE (hereafter 
SMILE-I), in order to study the background gamma 
rays for observations of celestial objects and verify the 
background rejection capability, we observed both atmo- 
spheric gamma rays, generated by the interaction be- 



tween cosmic-ray particles and nuclei in the atmosphere, 
and diffu se cosmic gam ma rays at an altitude of up to 
- 35 km (lTakadall2007h . 

Since the 1960s, several experiments have been con- 
ducted to observe diffuse cosmic gamma rays. The first 
detection was using th e Ranger 3 satellite be tween 70 keV 
and 1.2 MeV in 1962 (jMetziger et al.lll964[ ). Thereafter, 
the existence of these rays was confir med using the sate l- 
lites ERS-18hetween 0.25 and 6 M eV (iVette et al.|[T970l) 
OSO-3 between 7 and 100 keV ([Schwartz fc PetersonI 
I19741). Avollo 15 in the energy range of 0.3-27 MeV 
(|Trombka eraI|ll973D. and K osmos 461 between 28 keV 



and 4.1 MeV (iMazets et al.lll975D. and by usmg sev- 
eral balloons (iBleeker fc Deerenberd 119701 : lOamle et al.l 
[19711 IMakinol [T975t IKinzer et al.l I1978D . While these 
gamma-ray observatories used phoswich counters, ob- 
servations by using a balloon-borne double- Compton 
telescope i n the energy range above 1 MeV began in 
the 1970s l|Sch5nfelder et al.lll977l IT980l ). Subsequently, 
observations by using satellites have continued with, 
for exam ple, the High Energy Astronomy Obsery atory 
{HE AO; 'Mar shall et all figsot IKinzer et all IT991. th e 
Solar M aximum Mission (gM iwHWatanab e et al.lll999D 
CGR O (|Weidenspointner et al. 2000; Kappa dath et al.l 
[T99l) . and INTEGRAL (fChurazov et al.. '.20071) that is 
based on coded aperture imaging. An MeV bump, which 
is in excess near a few MeV in the Apollo observations, 
remained a mystery for a long time but is now rejected, 
because COMPTEL revealed that it was an artifact from 
delayed gamma -rays produced in radioactiv e decays in 
the instrument (jWeidenspointner et al.ll2000l ). 

Atmospheric gamma rays have been observed by bal- 
loon experiments s i nce the 1960s (iPeterson et al.lir972l: 
Kinzer et all 119741 : iRvan et al.l 119771 : iSchonfelder et al.l 
1980^ Low-energy atmospheric gamma rays are 



mainly produce d by bremsstrahlung from secondar y cas- 
cade electrons ([Peterson et al.l 119731 : iDaniol I1972D , and 
the electron-po s itron annih i lation line is also detected 
([Peterson! Il963t iLing et all I1977D . Some models for 
atmosp heric gamm a rays exist, for example, the Ling 
model (Lingl Il975l ) an d the Qine tiQ Atmospheric Radi- 
ation Model (QARM; iLei et al.l[200& ). Previous obser- 
vations showed that the fiux of atmospheric gamma rays 
depends on the zenith angle, and models produced by 
the following scientists are based on the observational 
results: Ling, Costa (Cost a et al.l Il984| ). and Graser 
(Graser & Schonfelder 1977). 

Cosmic and atmospheric charged particles and neu- 
trons could present large backgrounds in gamma-ray ob- 
servations; thus, th e estimation of the fluxes of these 
parti cles is essential ([Wunderer et al.l [20061 : iBowen et al.l 
l2007f l. The fluxes of these particles fluctuate on the basis 
of solar activity, observation altitude, and cutoff rigidity. 
A few models for energy spectra of these particles at 
balloon altitudes have been proposed, such as Mizuno's 
model ([Mizuno et al.ll2004[ ) and QARM. 

Recently, in addition to an ETCC, balloon experiments 
with new detection methods have been performed to 
achieve higher sensitivity in the sub-MeV/MeV band. 
For example, a double-Compton camera using a liq - 
uid Xe detector: the LXe GRIT (lAprile et all [2001: 
multiple Compton cameras ([Kamae et al.lll98'7l) by us- 
ing semiconductor s: th e nuclear Compton telescope 
fNCT: .Boggs ct al.,i2007f ). a Si/CdTe Compton telescope 
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(|Takahashill2005[ ) , and the Tracking and Ima ging Gamma 
Ray Experiment (TIGRE: [Z^ch et al.|[2008l) : an electron- 
tracking Compton camera using an electron tracker made 
of a semiconductor: the Medi um Energy Gamma -ray As- 
tronomy (MEGA) telescope fBloser et al."2006'); and a 
gamma-ray lens: CLAIRE fivon Ballmoo s ct al, 2005]) . 
Also, such experiments have renewed the attention given 
to observations of diffuse cosmic gamma rays, atmo- 
spheric gamma rays, and the background formed by cos- 
mic rays at balloon altitudes. 

In this paper, we report the observational results of 
SMILE-I. We first reconstructed the gamma-ray events 
and then estimated the background produced by the 
interaction between cosmic ray and the materials sur- 
rounding the detector during the flight. Next, we simu- 
lated the scattering of diffuse cosmic gamma rays in the 
atmosphere. After that we obtained the fluxes of dif- 
fuse cosmic and atmospheric gamma rays measured with 
SMILE-I. Finally, we determined the detection sensitiv- 
ity of SMILE-I, and estimated the expected detection 
sensitivity of a future SMILE mission. 




Fig. 1. — Schematic view of our electron-tracking Compton cam- 
era (ETCC). The ETCC consists of a gaseous tracker that detects 
the Compton-recoil electrons and a GSOiCe pixel scintillator ar- 
ray (PSA) that detects the Compton-scattered gamma rays. The 
distance between the drift plane and the gas electron multiplier 
(GEM) is 14 cm, and that between the GEM and the micro-pixel 
chamber (/i-PIC) is 0.5 cm. The effective area of the time projec- 
tion chamber (TPC) itself is the area between the drift plane and 
GEM. 



2. INSTRUMENT 

Figure [T] shows a schematic of the ETCC developed by 
us. The electron tracker detects the three-dimensional 
track and the energy of the Compton-recoil electron, 
while the absorber detects the absorption point and en- 
ergy of the Compton-scattered gamma ray. By summing 
the momenta of both the recoil electron and the scattered 
gamma ray, we obtain the momentum of the incident 
gamma ray as 

-Pi', (1) 



where Pq, pt^, and are the four-dimensional momenta 
of the incident gamma ray, scattered gamma ray, and 
recoil electron, respectively. Thus, we can obtain a com- 
pletely ray-traced gamma-ray image. The unit vector of 
the incident gamma ray r is described by 



sin (f) 
tan a 



g 



sin 4> , 



-g 



E^+K, ' 



(2) 
(3) 



where g and e are unit vectors in the directions of the 
scattered gamma ray and the recoil electron, respectively, 
and (j) is the scattering angle given by 



1 



rrieC 



1 



(4) 



where E^, Ke, We, and c are the energy of the scat- 
tered gamma ray, kinetic energy of the Compton-recoil 
electron, electron mass, and light speed, respectively, a, 
which is the angle between the scattering direction and 
the recoil direction, as shown in Fig. [TJ is geometrically 
measured for each incident gamma ray. 



cos Ckgeo = g • e. 



(5) 



and this angle is also obtained from a calculation by the 
Compton scattering kinematics: 



cos aki' 



1 



Ke + 2meC2 ■ 



(6) 



By comparing a^eo with a^- 



^ cos OLgeo 



cos aki 



(7) 



where Aq, is a cut parameter, we can select only Compton 
scattering events. Thus, Compton imaging with electron 
tracking is an effective method for MeV gamma-ray as- 
tronomy, with the serious problem of a large background, 
as described in Section 1. 

For SMILE-I, we have constructed an ETCC with a 
detection energy range 100 keV-1 MeV. Our ETCC con- 
sists of a gaseous tracker for detecting Compton-recoil 
electrons by using gas avalanche detectors to read it. In 
addition, it has a pixel scintillator array (PSA) for de- 
tecting Compton-scattered gamma rays. We provide a 
detailed description of our ETCC design in the following 
paragraphs. As a tracker for detecting Compton-recoil 
electrons, we developed a gaseou s time pro jectio n cham- 
ber (TPC; ,Kubo et al.i r2003: Miuchi et al.ll2003[ ) with a 
volume of 10x10x14 cm'^ filled with a gas mixture of 80% 
Xe, 18% Ar, and 2% C2H6 in mass ratio and sealed at 1 
atm. The readout of the TPC consists of g as avalanche 
detectors: a gas electr on multiplier (GEM; iSaulil 119971 : 



iTamagawa et al" 



PIC: lOchi et~ 



200' 



20011 



and a micro-pix el chamber (^- 
iNagavoshil [20041 ). The latter is 
our original gaseous two-dimensional imaging detector 
with micro-pixel electrodes produced using printed cir- 
cuit board technology. A seed electron drifts into the 
50-/xm-thick GEM with a velocity of 2.5 cm fis~-^ at an 
electric field of 400 V cm~^, and the first multiplication 
is caused by the GEM with a gain of approximately 10. 
Then, the multiplied electrons drift to the ^-PIC, and 
the second multiplication is caused by the /x-PIC with a 
gain of 3 x 10^. Therefore, we obtain a high gas gain 
above 3 x lO''^, which is enough to detect Compton-recoil 
electrons. The signals from the /i-P IC are read using 
ampl ifier-shaper-discriminator chips (jSasaki Sz Yoshidal 
Il999f) . which feed an output signal from a charge am- 
plifier to a flash analog-to-digital converter to measure 
the recoil-electron energy, and a discriminated digital 
signal to a position encoder that encodes the electron 
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Fig. 2. 
ETCC. 



Cross-sectional view of a mass model of the SMILE-I 



track by using field-programmab le gate arrays with a 
100-MHz clock (iKubo et al.ll2005D . Because the /z-PIC 
consists of fine-structure pixels with a pitch of 400 ^m, 
this TPC has a fine three-dimensional position resolu- 
tion of ,, / (T^[^ + (Tn = 500 /im, where a± and an are reso- 



lutions perpendicular and parallel to the drift direction, 
respectively (Takada 2007). As an absorber for detect- 
ing Compton-scat tered gamma rays, we used 33-pixel- 
scintillator arrays (jNishimura et al.ll2007f ). each of which 
consists of 8 X 8 GSO:Ce scintillator pixels with a pixel 
size of 6 X 6 X 13 mm^. For the photon sensor of the scin- 
tillation camera, we selected a multi-anode photomulti- 
plier tube (Hamamatsu Photonics, Flat-Panel H8500), 
which consists of 8 x 8 anode pixels with a pixel size of 
6x6 mm^ . With chained resistors connecting the anodes 
to reduce the number of readout channels, we obtained 
the position of a hit pixel by the charge-division method 
(jSekiva et al.ll2006t ). A plastic scintillator with a size of 
30 X 30 X 0.3 cm^ was placed 21 cm above the GEM as 
a veto counter for reducing triggers by charged particles. 
When an energy deposit in the GSO scintillators was over 
30 keV, a trigger was generated and the data-acquisition 
system waited for 8 /is (the maximum time for seed elec- 
trons to drift to the /i-PIC) for the TPC signals. When a 
trigger from the veto counter was generated, the system 
waited for 100 ^s (the recovery time of the undershoot 
of the large signal caused by charged particles) after the 
veto trigger. 

To obtain the fluxes of diffuse cosmic gamma rays and 
atmospheric gamma rays by SMILE-I, the detection effi- 
ciency and the effective area for parallel light is required. 
However, irradiating the entire system of SMILE-I by 
using parallel light in a laboratory is very difficult. In- 
stead, w e simulated the irradiat ion with Geant4 (ver 9.0- 
patchOl: lAgostinelh et alll2003[ ) by using a mass model 
shown in Fig. [5] In the ETCC simulator, the absorber 
had an energy resolution of 12 x (Energy/662keV)~°'^^ 
% at full width at half maximum (FWHM) and an en- 
ergy threshold of 30 kcV, the tracker had an energy res- 
olution of 45 X (Energy/32keV)"°-^^ % at FWHM and 
an energy threshold of ~ 3 keV, and the veto counter 
had an energy resolution of 50% at FWHM and an en- 
ergy threshold of 300 keV. The time resolutions of the 
absorber and the veto counter were not considered, be- 
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Fig. 3. — Dependence of detection efficiency on the energy of 
incident gamma rays for a point source near the detector, obtained 
by experiments (circles) and simulation (triangles). The depen- 
dence of efficiency on the incident energy for parallel light at a 
zero degree zenith angle from the Geant4 simulation using a mass 
model of Fig. [2] is also shown (solid line). The error bars represent 
only statistical uncertainties. 
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Fig. 4. — Effective area as a function of zenith angle for a point 
source near the detector obtained by experiments (filled circles) 
and simulations (open circles), and for parallel light with energies 
of 150 keV (dashed line), 600 keV (solid line), and 1500 keV (dotted 
line), obtained by Geant4 simulation using a mass model shown in 
Fig. [2] The error bars represent only statistical uncertainties. 



cause our ETCC did not use the event selection based on 
the time-of-fiight between the Compton scattering tar- 
get and the absorber. We obtained the detection effi- 
ciency of gamma rays reconstructed within zenith angles 
below 60 degrees for irradiation from the zenith direc- 
tion as a function of the incident energy, as shown in 
Fig. [21 the detection efficiency was approximately 10"'* 
between 150 keV and 1.5 MeV. To estimate the system- 
atic uncertainty of the detection efficiency, we measured 
and also simulated the detection efficiency for a radioiso- 
tope under the condition of gamma-ray irradiation from 
a point source placed at a distance of 50 cm from the 
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drift plane of the TPC. The difference between the ex- 
perimental measurement and simulated data was approx- 
imately 15%, as shown in Fig. |3l and we incorporated 
this difference into the systematic uncertainty for the 
detection efhciency. In order to verify our mass model 
simulation of parallel light and to determine if the detec- 
tion efficiency uncertainty depends on the zenith angle, 
we plotted the effective area obtained by simulation for 
parallel light between 150 keV and 1.5 MeV as a func- 
tion of zenith angle (Fig. |4]) . For the verification of the 
simulator, we measured and also simulated the effective 
area for a point radioisotope. The difference between 
the simulation and experimental data does not show any 
systematic dependence on incident direction. Thus, we 
verified our mass model simulation, and assumed that 
the systematic uncertainty of the detection efficiency is 
not dependent on the zenith angle. Figure |4] also shows 
that the acceptance field-of-view (FOV) of SMILE-I was 
120 degrees (3 sr) full width at half maximum (FWHM) 
for 150 keV, 200 degrees (7 sr) for 600 keV, and more 
than 7 sr for 1500 keV, respectively, which was wider 
than that of COMPTEL that was 2 sr (FWH M) at 1.2 
MeV and 0.5 sr at 6.1 MeV (Sch6nfcldereLalj[l993) . Ta- 
ble [T] shows the energy and the angular resolutions of the 
ETCC as a function of the incident energy, measured at a 
ground-based experiment. The angular resolution of the 
ETCC is defined by two parameters: the angular resolu- 
tion measure (ARM), which is the accuracy of the scat- 
tering angle 0, and the scatter plane deviation (SPD), 
which is the det ermination accura cy of the Compton- 
scattering plane ()Bloser et al.ll2002l ). In general, the the- 
oretical limit of the ARM is Doppler broadening caused 
by the undetectable momentum of the Compton-target 
electron, and the SPD is limited by the multiple scatter- 
ing of the Compton-recoil electron; thus, an ETCC using 
a gas with a higher atomic number has a poorer angular 
resolution due to the larger multiple scattering, although 
it has a higher detection efficiency. The ARM resolution 
due to Doppler broadening is 3.3 degrees for 200 keV and 
0.7 degrees for 1 MeV at FWHM. The SPD resolution 
due to multiple scattering of the recoil electron with re- 
coil energy of tens of keV is approximately 100 degrees 
at FWHM. In fact, an ETCC using an Ar gas had the 
ARM and SPD of 8.4 degrees and 89 degrees (FWHM), 
respectively, at 662 keV (jTakada 2007), and these values 
were better than those obtained using the higher atomic 
number Xe gas in the SMILE-I ETCC detector. Because 
we intended to observe not a celestial point source but 
diffuse cosmic gamma rays and atmospheric gamma rays 
in SMILE-I, we gave preference to the higher efficiency 
over the better angular resolution and, thus, adopted Xe 
gas for our ETCC. 

For the balloon flight, the ETCC was placed in an alu- 
minum vessel with a diameter of 1 m, a height of 1.4 m, 
and a thickness of 3 mm. The vessel was maintained at 
1 atm and the vessel was fixed to an aluminum gondola 
with a size of 1.2 x 1.5 x 1.6 m'^. On two sides of the 
gondola, batteries and ballast boxes were attached, and 
the gondola was packed with expanded polystyrene. In 
addition, the SMILE-I gondola had a fine pressure gauge 
for measurement of atmospheric pressure, a global posi- 
tioning system receiver for measurement of the balloon's 
altitude and geographic position, and two clinometers 
and two geomagnetic aspectmeters for determining atti- 



TABLE 1 

Energy resolution and angular resolutions 
(given by the FWHM) of the ETCC 



Energy 


Energy resolution 


ARM 


SPD 


[keV] 


[%] 


[degrees] 


[degrees] 


166 


37 ±2 


46 ±2 


189 ± 17 


356 


20 ±2 


24 ±2 


181 ± 10 


511 


15 ± 1 


21 ±2 


180 ± 15 


662 


14 ± 1 


18 ± 1 


183 ± 13 


835 


15 ± 2 


21 ±4 


185 ± 17 



tude. 

3. BALLOON FLIGHT 

The SMILE-I balloon was launched from the Sanriku 
Balloon Center of the Institute of Space and Astro- 
nautical Science/ Japan Aerospace Exploration Agency 
(ISAS/JAXA; 39.16°N, 141.82°E) on September 1, 2006 
at 06:11 Japan Standard Time (JST). The flight path 
is shown in Fig. [5j Figure |6] shows the time variations 
of the balloon's altitude and atmospheric pressure. At 
08:56, the balloon reached an altitude of 35.0 km and 
began level flight. From 08:56 to 10:15, the altitude was 
constant at 35 km and the atmospheric pressure was 5.4 
hPa (5.5 g cm~2)^ and from 11:20 to 13:00, the altitude 
was 32 km and the pressure was 8.5 hPa (8.7 g cm~^). 
During level flight between 08:56 and 13:00, the time 
average of the atmospheric depth was 7.0 g cm^^, the 
live time of the observation was 3.0 hr, and the cutoff 
rigidity calculated by QARM was 9.7 GV. Between 12:06 
and 12:33, we suspended gamma-ray observation, and 
the camera was operated in the charged-particle track- 
ing mode, which was used to check the performance of 
the tracker with a trigger of any two hits on GSO pixels. 
At 12:33, we resumed gamma-ray observation. Then, we 
turned off the system power at 12:59. Finally, we cut the 
gondola from the balloon at 13:20. The gondola landed 
in the sea at 13:45, and we successfully recovered it at 
14:32. During the whole flight, no serious problem oc- 
curred with the balloon system or the detector. 

4. RESULTS & DISCUSSION 

4.1. Event Reconstruction 

In this flight, 2.2 x 10^ events were acquired during 
the 6.8 hr after the launch. If an incident gamma ray 
had been scattered in the TPC and then hit more than 
two pixels in the GSO:Ce absorber, we could not re- 
construct the gamma-ray event because the sequence of 
interactions in the absorber was unknown. We, there- 
fore, selected events with a single hit in the absorber, 
and 1.1 X 10^ events remained. Figure [7] shows the en- 
ergy deposit in the TPC and the track length of those 
events. If a Compton-recoil electron stops in the TPC, 
the event is plotted around the dashed line in this flg- 
ure, which represents the relationship obtained from the 
Geant4 simulation. On the other hand, if the electron 
escaped from the TPC, the event is located around the 
solid line, which represents the relationship calculated 
from energy deposit per unit length for minimum ion- 
izing particles. This is because with increasing energy 
of the recoil electron its energy deposit per unit length 
of the electron escaping from the TPC approaches that 
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Fig. 5. — Flight path of SMILE-I. The balloon was launched 
from the Sanriku Balloon Center (SBC) of ISAS/JAXA in Japan 
(39.16N, 141. 82E), and it landed in the sea approximately 30 km 
from the SBC. 
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Fig. 7. — Track length and deposited energy in the TPC for 
events that occurred during level flight (scatter plot). The solid 
and dashed lines represent the relationship calculated from energy 
deposit per unit length of minimum ionizing particles and the rela- 
tionship simulated using Geant4 for electrons stopped in the TPC, 
respectively. The hatched area represents the event selection de- 
scribed by equations ([Sj and (|9]l. 




Fig. 6. — Time variation of the balloon's altitude and atmo- 
spheric pressure. 



of minimum ionizing particles. Because the energy of a 
Compton-recoil electron is necessary for the gamma-ray 
reconstruction, 6.5 x 10'^ events remained after we se- 
lected events in the case that the electron stopped in the 
TPC. We used the following criteria to select events in 
which the electron stopped in the TPC: 



- 1.8 X 10^3 



[kcV] 



<2, 



K,>15 [keVl 



(8) 
(9) 

where is the length of the detected track in the tracker. 
We know from our experience with ground-based exper- 
iments using the prototype trackers that an electron de- 
posits all its energy in the TPC if it satisfies equation ([5]). 
The energy spectrum of the TPC has a line component 
at 8 keV which corresponds to the characteristic CuKa 
X-ray line from the copper electrode of the TPC. Be- 



cause these events are not Compton events, we used the 
threshold given in equation (jH]) to reject them. Thus, we 
selected Compton events in which the electron stopped 
in the TPC in order to be able to reconstruct the incident 
gamma rays. Finally, we selected events whose Compton 
scattering occurred in the fiducial volume of 9 x 9 x 13 
cm^ by using the angle a and equation ([7|) with Aq, of 



A„ = 0.075 X 



Ke \Eg - mgC^I 



(10) 



2.1 X 10"^ events then remained. This criteria of the angle 
a is equivalent to |1— |rp| < 0.15, where r was described 
by equation ([3]) and |r| is equal to unity under the con- 
dition of cos Q!geo — cos akin ■ The time variation of the 
count rate of the reconstructed events is shown in Fig. [8l 
which has a maximum near the Pfotzer maximum and 
is thereafter constant during level flight, with a stan- 
dard deviation of 20%. During level flight, we obtained 
881 reconstructed events in the detected energy range of 
100 keV-1 MeV, 420 events of which were downward at 
zenith angle less than 60 degrees. The energy spectra of 
both the downward events and the events in all directions 
are shown in Fig. [S] 

4.2. Instrumental Background Simulation 

In a balloon-borne experiment, particles incident on 
a detector are not only gamma rays but also protons, 
electrons, neutrons, and others. These particles and the 
gamma rays induced by interactions between these parti- 
cles and the pressurized vessel or materials surrounding 
the detector could form a part of the background. In 
particular, a neutron causes elastic scattering, so that 
when a neutron interacts with both the scatterer and 
the absorber of a Compton camera, the neutron event 
may be incorrectly reconstructed as a Compton-scattered 
gamma ray in Compton- type detectors. The electron 
tracker of an ETCC can distinguish between electrons 
stopped in the TPC and other charged particles because 
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Fig. 8. — Time variation of the rate of reconstruction of events. 
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Fig. 9. — Energy spectra of the events that were reconstructed 
as gamma rays during level flight in all directions (circles), and at 
zenith angles below 60 degrees (squares). 
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Fig. 10. — Mass model for the instrumental background simula- 
tion. The model of the electron-tracking Compton camera (ETCC) 
is the same as in Fig. [2] 



of the difference in the track length as a function of 
the kinetic energy, as shown in Fig. [T] The purpose of 
the SMILE-I experiment was the measurement of dif- 
fuse cosmic and atmospheric gamma rays. Therefore, 
atmospheric gamma rays produced in the atmosphere by 
cosmic-ray interactions are considered not as background 
but as signals in this experiment. 

To estimate background events, we calculated the 
background radiation at the balloon's altitude by using 



a Geant4 simulation and a mass model shown in Fig. 1101 
The balloon gondola was not included in the mass model, 
because most of the parts of the balloon gondola were 
expanded polystyrene and the effect on secondary par- 
ticle production was negligible. Due to the low density 
of the electronics boxes and the large distance between 
the ETCC and the battery, the effect of the instrumen- 
tal background can be neglected. Therefore, this mass 
model does not include a gondola, electronics box, and 
battery. The Geant4 hadronic processes were elastic and 
inelastic scattering for the low- and high-energy models, 
and absorption for pions, kaons, protons, anti-protons, 
and anti-neutrons; elastic and inelastic scattering, and 
capture with a high-precision model below 19 MeV for 
neutrons; elastic and inelastic scattering with a low- 
energy model for deutrons, tritons, and alpha particles. 
In addition, G^Decay and G4RadioactiveDecay processes 
were included for prompt processes and delayed photons 
from radioactive decays. We did not include the chance 
coincidences between triggers from different components 
of the radiation background, since their contribution is 
expected to be small by the background rejection method 
described in section 2. 

As described in section 1, there are several radiation 
models for the fluxes of incident particles at balloon al- 
titudes. Mizuno's model is a radiation model based on 
the previous observations of protons and alpha particles 
above 100 MeV, electrons and positrons above 100 MeV, 
gamma rays between 1 MeV and 100 GeV, and muons 
between 300 MeV and 20 GeV. However, below these en- 
ergy ranges, the flux of each of these particles is simply 
assumed to be a power-law spectrum with an index of 
— 1, and the neutron flux is not described at all in this 
model. In contrast, the QARM model, based on pre- 
vious observations and some simulations, has fluxes of 
neutrons in addition to protons, electrons, gamma rays, 
pions, and muons in wider energy ranges between 100 
keV and 100 GeV and includes the effect of solar mod- 
ulation on the spectral shape and normalization. The 
parameters of the QARM model are date, geographic 
position, altitude/residual pressure, magnetic field dis- 
turbance index, and the primary-particle spectrum. Be- 
cause the model of low-energy charged particles and neu- 
trons was important in sub-MeV gamma-ray observa- 
tions, we adopted the QARM spectra shown in Fig. [TTl 
assu ming a primary-particle spectrum of galactic cosmic 
rays (iLei et al, 2004,) und er a magnetic field disturbance 
index Kj, (|Mavaudlil980l ) of 3 and an atmospheric depth 
of 7 g cm~^. The zenith angle distribution of QARM for 
public use is divided in only two regions: the zenith an- 
gle range of 0-90 degrees, and 90-180 degrees. Figure [Til 
represents the downward fluxes in the zenith angle range 
of 0-90 degrees and the upward fluxes in the zenith angle 
range of 90-180 degrees. We estimated an uncertainty in 
the particle fluxes of QARM to be 20% by varying the 
atmospheric depth between 5.3 g cm~^ and 8.5 g cm~^ 
and the magnetic field disturbance index between 2 and 
4. 

For the simulation of the instrumental background 
coming from any direction, we assume the incident di- 
rection of the primary particles to be downward (zenith 
angle = degrees) or upward (zenith angle = 180 de- 
grees), because we do not have any observational data 



Fig. 11. — Energy spectra of protons (solid lines), electrons 
(dashed lines), and neutrons (dotted lines) for a background sim- 
ulation at an atmospheric depth of 7.0 g cm~^. Red shows the 
downward flux in the zenith angle range of 0-90 degrees and blue 
shows the upward flux in the zenith angle range of 90-180 degrees. 
All fluxes were calculated by QARM. 



Fig. 12. — Simulated source function of the secondary photons 
produced in the instrument by protons (dotted line), electrons 
(dashed line), and neutrons (long dashed line) at zenith an gles 
between and 90 degrees, spectra of which are shown in Fig. 1111 
at an atmospheric depth of 7.0 g cm~^. The solid line shows the 
total flux. 



of the horizontal fluxes or the zenith angle distributions 
of protons, neutrons, and electrons. For the incident 
fluxes, we adopted the QARM flux in the zenith angle 
range of 0-90 degrees (Fig. [TT|) as the downward pri- 
mary particle flux, and the QARM flux in the zenith an- 
gle range of 90-180 degrees (Fig. [TT]) as the upward pri- 
mary particle flux. The simulated spectra of instrumen- 
tal background gamma rays produced in the materials 
surrounding the detector are shown in Fig. 1121 The spec- 
tra shown in Fig. [12] are the flux averaged over the solid 
angle in the zenith angle range between and 90 degrees. 
Gamma rays induced by electrons dominate in the energy 
range of SMILE-I, and neutron-induced gamma rays con- 
tribute slightly at several MeV. In addition, there is a 
line component of electron-positron annihilation at 511 
keV. In order to estimate the contribution of secondary 
gamma rays, neutrons, and charged particles incoming 
to the gaseous tracker from all directions, we simulated 
their fluxes (Fig. [T3|). Throughout this paper, secondary 
gamma rays refer to photons due to the instrumental 
background only. Secondary gamma rays contribute the 
majority of the background, with the rest coming from 
neutrons, electrons, positrons, protons, and charged pi- 
ous. The rejection inefflciencies for background particles 
after Compton reconstruction under the constraints of 
equations ([5]), (O, and ([TU]) were estimated by the simu- 
lation to be 1.8 ±0.3 X 10"^ for neutrons, 3.4 ±0.2 x IQ-* 
for electrons, and 2.5±0.2xl0^'*for protons. Therefore, 
the simulation predicted that the background event rates 
of gamma rays, neutrons, and charged particles during 
level flight would be 7.4±2.6 x lO'^ s^^, 6.7±2.6 x 10"* 
s~\ and 8.7 ± 2.3 x 10"^ s"\ respectively. The un- 
certainties are caused by the uncertainty in the back- 
ground model. Because the rate of reconstructed events 
during the level flight of SMILE-I was 3.9 x 10"^ s'^ 
the contributions of background gamma rays, neutrons, 
and charged particles to the obtained events were 19%, 
1.7%, and 0.23%, respectively. Thus, we estimated that 
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Fig. 13. — Simulated differential fluxes of particles incident on 
the gaseous tracker from all directions for secondary gamma rays 
(open circles), neutrons (crosses), electrons (squares), positrons 
(asterisks), protons (filled circles), 7r+ (stars), and 7r~ (triangles), 
at an atmospheric dept h of 7.0 g cm~^. These spectra include both 
primary (shown in Fig. Illll and secondary particles. 

98.0 ± 0.7% of the reconstructed events were gamma- 
ray events, and the contribution of the other particles 
was 2.0 ± 0.7%, thanks to the powerful background re- 
jection of the ETCC. Figure[T4l shows the source function 
of instrumental secondary gamma rays incoming to the 
detector from the zenith direction as a function of atmo- 
spheric depth. These spectra do not include any event 
selection described in section 4.1. In later analysis, we 
use this flux of gamma rays as the background. 

4.3. Gamma-ray Scattering in the Atmosphere 

In observing diffuse cosmic gamma rays at a few hun- 
dred keV, there is a contribution of cosmic gamma rays 
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Fig. 14. — Differential flux of instrumental background photons 
produced by protons (dotted line), electrons (dashed lino), and 
neutrons (long dashed line) in the 125-1250 keV region, as a func- 
tion of atmospheric depth. These fluxes represent the background 
source function averaged in the zenith angle between degree and 
90 degrees at the sensitive area of the ETCC. The solid line showrs 
the total flux. 

that are scattered into the aperture of the detector af- 
ter one or more Compton scatterings in the overlying 
atmosphere. The cosmic gamma-ray flux jc{z, 9) at the 
atmospheric depth of z and the zenith angle of 9 can be 
separated into two components, 

Mz,e)^Uz,9) + Mz,e). (11) 

fn{z, 9) is the flux of cosmic gamma rays that enter the 
aperture of the detector without any interaction in the 
overlying atmosphere. In the zenith direction {9 = 0), 
the depth dependence of /„(z, 9) is given by 

fniz,0) = /„(0,0)exp(-Tfotz) = r]e:>q>{-Ttotz), (12) 

where 77 is the flux of cosmic gamma rays at the top of 
the atmosphere, and Ttot is the cross section of the total 
attenuation per unit mass. On the other hand, fsiz, 9) is 
the flux of cosmic gamma rays scattered into the aperture 
of the detector after one or more Compton scatterings in 
the overlying atmosphere. 

To obtain the ratio of the scattered component to the 
total diffuse cosmic gamma-ray flux, we simulated the 
transport of gamma rays in the atmosphere with Geant4. 
We defined 80 atmospheric layers from sea level to an alti- 
tude of 80 km, where each layer was a 1-km thick sphere, 
filled with air with the density appropriate for each al- 
titude. Initial gamma rays were generated at random 
positions on the surface of a sphere of 650 km radius and 
given random initial directions, and we traced their his- 
tories. Then, the ratio of the scattered component was 
calculated as a function of energy, atmospheric depth, 
and zenith angle. Figures [TSl and [TBI show the scattered 
component ratio A(z, 9) as a function of the detected en- 
ergy, 

assuming a photon index of 2.0 in the primary spec- 
trum. There are two breaks in Fig. 1151 The first 
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Fig. 15. — Ratio of the scattered component A of diffuse cos- 
mic gamma rays to the total (scattered plus nonscattered) gamma 
rays as a function of the energy at zenith angles between and 20 
degrees at an atmospheric depth of 8.0 g cm~^. The solid and dot- 
ted lines show our simulation results with and without Rayleigh 
scattering, respectively. The long-dashed line, dashed line, 
andj^tars show the simula tion bvporstman Hi Horstman-Morettil 
( 1971' ), the calculation by IMakind 119701) . and the simulation by 
Scho nfelder et aLl II1977I ). respectively. Each model is at an atmo- 
spheric depth of 8.0 g cm"'^. 




Fig. 16. — Our simulation results with Rayleigh scattering of 
the ratio of the scattered component A of diffuse cosmic gamma 
rays as a function of the energy at zenith angles between and 20 
degrees at atmospheric depths of 52.6, 11.5, 2.78, and 0.78 g cm~^ 
(upper to lower lines). There are five data sets for each curve: 10— 
100 keV, 100 keV-1 MeV, 500 keV-5 MeV, 1-10 MeV, and 2.5-25 
MeV, and the points in overlapped energy range are overlaid. 

break is at around 30 keV, where absorption and scat- 
tering dominate below and above the peak energy, re- 
spectively. The energy of the first break is constant at 
any atmospheric depth, as shown in Fig. [161 The sec- 
ond break in Fig. [15] appears at approximately 1 MeV. 
The energy of the second break decreases as the atmo- 
spheric depth decreases, as shown in Fig. [1^1 and the 
optical depth at the second break is approximately one. 
Therefore, for gamma rays with energies above the sec- 
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Fig. 17. — Scattering component ratio A for the initial spectra 
with power indices of 2.0 (solid line) and 2.5 (dashed line) at an 
atmospheric depth of 7.0 g cm~^ and zenith angles between and 
60 degrees. 



Fig. 18. — Ratio of the scattering component A as a function of 
the atmospheric depth at zenith angles between and 5 degrees, in 
the energy ranges 30-50 keV, 100-150 keV, 300-500 keV, 750-1000 
keV, 1000-1500 keV, 3000-5000 keV, and 5000-10000 keV. 



ond break, the probability of scattering is so low that 
A decreases. Moreover, Rayleigh scattering contributes 
a large amount (more than 10%) to the scattering com- 
ponent ratio under 50 keV. Compared with the previous 
calculation and simulations of the scattered component 
ratio function, our simulation result has a difference of 
<: 10% from th e simulation of Horstman at 30-100 keV 
(|Horstman fc H orstman-Moretti 197l|), ^ 20% from the 
calculation of Makino at hundreds of keV fM akinolHQTOl ). 
and ~ 10% from the simulat ion of Schonfelder above 1 
MeV (jSchonfelder et al.|[l977l ). Makino's calculation as- 
sumed that the scattering angle dependence of the cross 
section was negligible, and only one scattering occurred. 
Therefore, Makino's calculation differs considerably from 
our result in the energy range under 100 keV because 
of the scattering number limitation, and it also differs 
above 100 keV because it neglects the scattering angle 
dependence of the cross section. The contribution of cos- 
mic gamma rays that are scattered in the atmosphere is 
not negligible in the sub-MeV/MeV region; thus, we also 
used Geant4 to simulate the scattering component ra- 
tio as a function of energy (for various power indices), 
atmospheric depth, and zenith angle (Fig. [T71- [TO|) . Fig- 
ure [T7] shows that the ratios for different power indices 
differ by less than 5% in the range between tens of keV 
and a few MeV; thus, the choice of the power index does 
not strongly affect our results. Figure [TBI shows that the 
atmospheric depth dependence of A is similar to that 
of Makino's calculation, oc log(l -I- Ttotz), although 

there is a 10%-20% difference in the energy range of 100- 
500 keV. Figure [19] shows that although the zenith de- 
pendence of A below a few MeV is weak in the region 
of zenith angles less than 50 degrees, it changes abruptly 
above 70 degrees. Thus, the zenith angle dependence be- 
low a few MeV is not consistent with 1/ cos 6*, where 9 
is the zenith angle of the incident gamma ray. For the 
higher energies, A is similar to l/cos0. 

4.4. Diffuse Cosmic and Atmospheric Gamma-ray 
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Fig. 19. — Ratio of the scattering component A as a function of 
the zenith angle in the energy ranges of 30-50 keV, 100-150 keV, 
300-500 keV, 750-1000 keV, 1000-1500 keV, 3000-5000 keV, and 
5000—10000 keV. The dashed line is proportional to 1/ cos0. 

To study the origin of the reconstructed events, we di- 
vided the obtained data at atmospheric depths between 
5.2 g cm~^ and 168 g cm~^ into 10 segments as Table [2j 
Then, we obtained the gamma-ray flux from the zenith 
direction (Q = 0-10 degrees) for each time segment by 
analyzing the reconstructed events with a simulated re- 
sponse matrix for incident energies of 100, 150, 200, 300, 
400, 511, 600, 1000, and 1500 keV and incident zenith 
angles of 0, 20, 40, 60, 80, 120, 150, and 180 degrees. 
In the calculation, the zenith angle dependences of both 
diffuse cosmic and atmospheric gamma-ray fluxes must 
be considered, because our ETCC had a wide FOV of 
3 sr. For diffuse cosmic gamma rays, we simulated the 
zenith angle dependence as described in section 4.3. For 
the zenith angle distribution of combined cosmic and at- 
mospheric gamma rays, only three models exist and the 
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conditions under which they apply are hmited, as foUows. 

1) In Ling's model, the source function, which is the pro- 
duction rate of atmospheric gamma rays at the gener- 
ated position, is defined on the basis of observed data 
but has no dependence on the zenith angle. In trans- 
port, only absorption is included, and scattering is not 
considered. Ling's model also describes the flux of dif- 
fuse cosmic gamma rays between 300 keV and 10 MeV 
at atmospheric depths smaller than 500 g cm~^, which 
is comparable to the results of several experiments, al- 
though the upward gamma-ray flux differs by a factor of 
three from the observed data (jSchonfelder et al.l 119771 ). 

2) Graser's model calculates atmospheric gamma rays 
based on bremsstrahlung using some models of the elec- 
tron flux. This model ignores scattering in transport 
and Tr'^-decay gamma rays. Thus, this model is effective 
in the energy range of 1-30 MeV. 3) Costa's model is an 
empirical model for the sum of cosmi c and atmospheric 
gamnra rays bas ed on observa.tions by iSchonfelder et al.l 
(|l977l[T980l ) and lRvan et"all (flQTl . This model fits the 
experimental data well but is applicable only at atmo- 
spheric depths under 10 g cm~^. Given the limitations 
of each of these models, it is difficult to choose which one 
best describes the zenith angle distribution of combined 
cosmic and atmospheric gamma rays. 

As described in section 4.3, the contribution of the 
scattering of cosmic gamma rays in the atmosphere is 
not negligible within the energy range of SMILE-I. Sim- 
ilarly, the scattering of atmospheric gamma rays in the 
atmosphere is not negligible. However, Graser's model 
ignores scattering in the atmosphere. Ling's model also 
does not include any scattering in transport. However, 
the difference from the observed zenith-angle dependence 
is less than that of Graser's model, because the source 
function of Ling's model was fitted to the sum fluxes 
of cosmic and atmospheric gamma rays obtained from 
previous observations. Moreover, Ling's model covers 
the range of atmospheric depths corresponding to all al- 
titudes of the SMILE-I flight. We, therefore, adopted 
Ling's model including the cosmic component. However, 
the Ling's model describes the flux only above 300 keV. 
Thus, we extended Ling's model, which was extrapolated 
to less than 300 keV with a single power-law spectrum flt- 
ted between 300 keV and 5 MeV. Below an atmospheric 
depth of 10 g cm~^, the difference between the extended 
Ling's model and Costa's model in the energy range be- 
low 300 keV is less than ±5% at all zenith angles. How- 
ever, we could not flnd a report on any observational 
data below 300 keV that could be compared with the 
extended Ling's model above the atmospheric depth of 
10 g cm~^. Therefore, although we adopted an extended 
Ling's model, we could not estimate its uncertainty. 

We used this extended Ling's model of the cosmic and 
atmospheric gamma rays along with the simulated re- 
sponse matrices to unfold the observed events with re- 
constructed directions in the zenith angle below 60 de- 
grees, resulting in the gamma-ray fluxes listed in Ta- 
ble [21 In this analysis, we assumed that the zenith angle 
dependence of the sum of diffuse cosmic, atmospheric, 
and instrumental gamma rays was the same as the ex- 
tended Ling's model. Using this data, the growth curve, 
which is the dependence of the gamma-ray flux on atmo- 
spheric depth, in the energy range of 125-1250 keV was 
obtained, as shown in Fig. 1201 For comparison with the 
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Fig. 20. — Reconstructed gamma-ray flux in the zenith angle 
range of to 10 degrees in the energy range 125-1250 keV as a 
function of atmospheric depth. The circles represent data taken 
with the balloon ascending. The dashed Une sho ws the sum of 
1) the simulated background (solid hne in Fig. I14|l . 2) the diffuse 
cosmic gamma-ray flux of our extended Ling's model between 125 
and 1250 keV, and 3) the atmospheric gamma-ray flux of the same 
model. The atmospheric component of the extended Ling's model 
is scaled by the cutoff rigidity. 

extended Ling's model, the atmospheric gamma-ray flux 
was considered to have a dependence on _R ~,V^'^, where 
Rcut is the cutoff rigidity ([Thompson fc SimpsonlllQSil ). 
Because the cutoff rigidities of the SMILE-I flight and 
Ling's model are 9.7 GV and 4.5 GV, respectively, the 
sum of the atmospheric gamma-ray flux of the extended 
Ling's model multiplied by (9.7/4.5)~^'^^ and the diffuse 
cosmic gamma-ray flux of the extended Ling's model be- 
tween 125 keV and 1250 keV are also shown in Fig. [20l 
The SMILE-I result was consistent with the sum of the 
extended Ling's atmospheric and cosmic model and the 
background gamma-ray flux represented by the solid line 
in Fig.[T4l At atmospheric depths from 5 to 100 g cm~^, 
the variability of the ratio of the background particle 
flux to the cosmic and atmospheric gamma-ray flux is 
slight (< ±5%). Thus, the ratio of the non-Compton- 
scattering events of neutrons or charged particles to the 
reconstructed events is approximately equal to that of 
the level flight, 2%, which was estimated in section 4.2. 
In other words, 98% of the reconstructed events were 
induced by cosmic, atmospheric, and instrumental back- 
ground gamma rays. Therefore, we calculated the fluxes 
of cosmic and atmospheric gamma rays by using the 
fluxes obtained from the SMILE-I observations assum- 
ing that the reconstruction of non Compton-scattering 
events was negligible. 

To obtain the fluxes of diffuse cosmic gamma rays and 
atmospheric gamma rays, we used the difference in the 
atmospheric depth dependence of each component. For 
the cosmic gamma-ray flux at an atmospheric depth of 
z, we assumed 



1 - A(z,0) 



exp(-Ttotz), 



(14) 



with equations ([TT]), ([T2|), and at = 0. 

On the other hand, the atmospheric gamma-ray flux 
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TABLE 2 

Obtained gamma-ray fluxes in the zenith angle range of 0-10 degrees for each time 

SEGMENT. 





Time 




Atmospheric depth 


Flux 


10 * photons 


— 1 —9 — 1 

s cm sr 


keV-i] 




[JST] 




[g cm-2] 


125-250 keV 


250-550 keV 


550-1250 keV 


125-1250 keV 


07 


20-07 


47 


145 (121-168) 


8.8 ± 1.2 


5.5 ± 0.6 


1.6 ±0.2 


3.4 ±0.3 


07 


47-08 


00 


95 (70-121) 


18 ± 10 


9.0 ± 1.2 


2.6 ±0.4 


6.0 ± 1.2 


08 


00-08 


12 


55 (40-70) 


12 ±2 


7.4 ± 1.0 


2.1 ±0.4 


4.6 ±0.5 


08 


12-08 


24 


32 (23-40) 


7.1 ± 1.1 


4.5 ± 0.5 


1.3 ±0.2 


2.8 ±0.3 


08 


24-08 


56 


12 (5.3-24) 


5.1 ±0.7 


3.2 ± 0.4 


0.92 ±0.13 


2.0 ±0.2 


08 


56-10 


08 


5.4 (5.2-5.6) 


3.0 ±0.4 


1.9 ±0.2 


0.53 ±0.07 


1.1 ±0.1 


10 


08-11 


00 


6.8 (5.4-7.9) 


3.9 ±0.5 


2.4 ±0.3 


0.69 ±0.09 


1.5±0.1 


11 


00-11 


30 


8.2 (7.7-8.6) 


3.7 ±0.5 


2.3 ±0.3 


0.65 ±0.10 


1.4 ±0.1 


11 


30-12 


06 


8.5 (8.3-8.7) 


3.8 ±0.6 


2.4 ±0.3 


0.67 ±0.10 


1.5±0.1 


12 


33-12 


59 


8.3 (8.0-8.7) 


3.2 ±2.3 


2.0 ± 1.2 


0.57 ±0.33 


1.2 ±0.5 



TABLE 3 

Obtained diffuse cosmic and atmospheric gamma-ray fluxes in the zenith 
direction (0-10 degrees). 



Energy 
[keV] 


Cosmic gamma rays^ 
[ph cm"^ sr"^ kcV^^] 


Atmospheric gamma rayst 
[ph cm-^ sr--*- keV-'" (g cm-^)--'] 


125-250 
250-550 
550-1250 


2.1 ±0.7tJ;° X 10-* 
8.3 ± 3.4t^ | X 10-5 
3.3±1.2l^;^ X 10-5 


2.0 ± o.etoi X 10-5 

7.8 ± 2.915 5 X 10-*^ 
3.3± l.Ot?;^ X 10-*^ 



t The first uncertainty term is statistical and the second is systematic. 



fajz) is nearly proport ional to the atmospheric depth 
(|Schonfelder et al.l[l977l) . 



Then, we fitted the observed growth curve to 



(15) 



1 - \{z,0) 



eiip{-Ttotz)+fb{z)+F5ii{z)S{EQ- 

(16) 

where ft is the flux of background photons estimated by 
the Geant4 simulation fFig.rH]). F511 is the flux of the at- 
mospheric annihi lation line described by equation (1) in 
ILing et all (|1977f ) , and C and 77' are the attenuated fluxes 
through the pressurized vessel. The attenuated fluxes 
are given as ^' = ^ exp{—Tabsdl) and 77' = 'r]exp{—Tabsdl), 
where Tabs is the cross section of photoelectric absorp- 
tion in aluminum, d = 2.7 g cm-^ and I = 3 mm are 
the density and the thickness of the pressurized vessel, 
respectively. The growth curves with the fitting func- 
tions are shown in Fig. dTJ and these functions serve to 
isolate individual contributions to the total flux so that 
separate fluxes for diffuse cosmic gamma rays and atmo- 
spheric gamma rays can be obtained. However, at larger 
atmospheric depths than approximately 100 g cm-^, we 
have a problem because the gamma-ray flux is not pro- 
portional to the atmospheric depth, as shown by Ling's 
model in Fig. [201 Therefore, the data at the highest 
atmospheric depth was not used in the fitting. The ob- 
tained fiuxes of cosmic and atmospheric gamma rays are 
listed in Table [3] The systematic uncertainties in the 
cosmic and atmospheric gamma-ray fluxes are caused by 
the following: 1) The uncertainty of 15% in the detection 
efficiency (section 2). 2) The difference of 5% in the rel- 
ative ratio of the zenith angle distributions between the 
extended Ling's model and Costa's model (section 4.4). 



3) The uncertainty of 20% in the background gamma- 
ray flux obtained by varying the parameters in QARM 
(section 4.2). 4) The uncertainty of 5% in the correction 
factor of scattering in the atmosphere (section 4.3). Fig- 
ures [51] and [22] show the differential flux of diffuse cosmic 
gamma rays and atmospheric gamma rays, respectively, 
wit h re sults of previous observations for comparison. In 
511 felYHl fluxes fr om observations for which the FOV was 
not described in the references were not plotted. In the 
energy range of 125-250 keV, the flux of atmospheric 
gamma rays obtained by SMIL E-I had a difference o f 
3 sigma from the observation of [Peterson et al] (|1972D . 
However, other results for the diffuse cosmic and atmo- 
spheric fluxes obtained by SMILE-I were consistent with 
previous observations and the Ling and QARM models 
of atmospheric gamma rays within 2 sigma. 

4.5. Sensitivity and Future Prospects 

Using the SMILE-I result, we estimated the contin- 
uum detection sensitivity of this camera. Diffuse cos- 
mic and atmospheric gamma rays, and those induced 
by charged particles and neutrons, which were observed 
in the SMILE-I experiment, would form the background 
when we observe a celestial point source. The minimal 
detectable flux Fmin at the significance of 3a is described 
as 



F„ 



If AE An 

^ef f Tabs 



(17) 



where / is the energy spectrum of background radiation 
in units of photons s~^ cm-^ sr-^ MeV~^, and A^ff, 
Tabs, AE, and Afl are the effective area, effective oId- 
servation time, en ergy resolution, a nd angular resolu- 
tion, respectively (iSchonfeldeii [200 ll) . From the results 
of the SMILE-I flight, the sensitivity of SMILE-I is ex- 
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pected to reach 4.1 X 10 ^ MeV s ^ cm ^ at 300 keV and 
1.2 X 10^2 jvjgv s-i cm-2 at 600 keV, assuming AE ^ E 
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Fig. 21. — Separation of atmospheric and cosmic components 
of gamma rays by fitting tfic growth curves in the energy ranges 
(a) 125-250 keV, (b) 250-550 keV, and (c) 550-1250 keV, in the 
zenith angle range of 0-10 degrees. The long-dashed, dotted, dot- 
dashed, and dashed lines represent the contributions from cosmic 
gamma rays, atmospheric gamma rays, instrumental background 
components, and the electron— positron annihilation line, respec- 
tively. The solid line represents their sum. 



and Tobs = 10" s (Fig. M 

For the next flight (hereafter SMILE-II) to observe the 
Crab Nebula from balloon altitudes, the ETCC needs 
upgrading in both the angular resolution and effective 
area. As described in equation the ARM, which is 
the accuracy of the scattering angle, depends mostly on 
the energy resolution of the absorber. Therefore, im- 
proving the energy resolution of the scintillation camera 
would give the ETCC better angular resolution. A good 
ARM resolution of 4 degrees (FWHM) at 662 keV was re- 
alized using a recently developed scintillator, LaBraiCe, 
with an excellent energy resolution of 3.0% and 5.8% at 
662 keV for a monolithic cry stal and an array of pixels, 
respectively (jKurosawa et al . 2009, 2010). From this re- 
sult, the expected ARM resolution would be 1.7 degrees 
and 1.3 degrees (FWHM) at 2 MeV and above 5 MeV, 
respectively. As for the effective area, we have devel- 
oped a larger ETC C based on a TPC with a volume of 
30 X 30 X 30 cm^ (jMiuchi et al.Jl2007D. and its g amma- 
ray imaging has been successful ( Ueno et aLll2008[ ) . With 
these improvements and the assumption of the instru- 
mental background scaled by the geometrical area of 
SMILE-I, it is expected that the ETCC for SMILE-II 
on a mission to detect the Crab Nebula during a level 
flight of 10"* s will have approximately ten times higher 
sensitivity than that of SMILE-I. 

We have additional ideas for improving the future 
ETCC. First, we are considering inserting an electron ab- 
sorber between the electron tracker and the scintillators. 
In the MeV region, because the energy of Compton-recoil 
electrons is as high as an MeV, the Compton-recoil elec- 
trons partly escape from the effective volume of the TPC. 
If an absorber that measures the energy of the escaping 
Compton-recoil electrons is placed between the TPC and 
scintillators, the sensitivity in the MeV region can be im- 
proved. Inserting the electron absorber has the added 
benefit of reducing the coincidence-timing window, be- 
cause we can take the coincidence between the electron 
absorber and scintillators. Second, we will adopt CF4 
gas as the electron tracker. Compared with Xe gas, the 
electron diffusion during electron drift and the multiple 
scattering effect of the recoil electron are both smaller. 
Thus, we will be able to trace the electron tracks more 
accurately. In addition, the position resolution of the 
Compton points and the accuracy of the recoil direction 
should be improved by updating the tracking algorithm. 
Last, we will develop the reconstruction algorithm on the 
basis of a reflexive calculation, such as the likelihood or 
maximum entropy method. If we realize an ETCC with 
a high ARM due to the use of LaBra PSAs with a pixel 
size of 6 X 6 X 40 mm"^, equipped with a TPC with a larger 
volume of 50 x 50 x 50 cm'' filled with CF4 gas at a higher 
pressure of 2 atm, the sensitivity will become ten times 
better than that of COMPTEL in the sub-MeV region, as 
shown in Fig. 1241 Moreover, when a photon with energy 
above 1.02 MeV causes electron-positron pair creation in 
the TPC, we can obtain the incident photon's momen- 
tum by measuring the m omenta of both the electron and 
positron with the TPC (jOrito et all 120031 ). Therefore, 
the camera will have better sensitivity above a few MeV 
than that shown in Fig. [24l In the near future, we will 
experiment on gamma-ray reconstruction by electron- 
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Fig. 22. — Differential fluxes of diffuse cosmic gamma rays as a function of energy. Open squares denote our results, and the uncertainty 
bars represent the la combined statistical and systematic uncertainties. References to previous observations are provided in the text. 



positron pair creation by using an ETCC and verify the 
result with a simulation, and also estimate the sensitivity 
of an ETCC onboard a satellite where the radiation envi- 
ronment is much different from that at balloon altitudes, 
by a performing the full-scale Monte Carlo simulation. 

5. CONCLUSION 

We are developing an ETCC as an MeV gamma-ray 
telescope to conduct an all-sky survey in the next genera- 
tion with sensitivity that is one order of magnitude higher 
than that of COMPTEL. To prepare for future observa- 
tions onboard a satellite, we performed a balloon exper- 
iment, SMILE. As the first step of SMILE (SMILE-I), 



we launched an ETCC consisting of a gaseous TPC with 
a volume of 10 x 10 x 14 cm'^ and GSO:Ce scintillation 
cameras onboard a balloon in 2006 for the observation 
of diffuse cosmic gamma rays and atmospheric gamma 
rays. In this flight, we successfully detected 2.1 x 10'^ 
reconstructed events from all directions, 420 of which 
were detected in an FOV of 3 sr during the live time 
of 3.0 hr of level flight. To evaluate the instrumental 
background, we estimated particles incident on the de- 
tector with simulations using QARM and Geant4; the re- 
constructed events included 1.7% neutron events, 0.23% 
charged-particle events, and 98% gamma-ray events. Be- 
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Fig. 23. — Differential fluxes of downward atmospheric gamma rays in the zenith direction (0—10 degrees), multiplied by {i?cut/9.7 GV)^'^'' 
as a function of energy. Open squares denote our results, and the uncertainty bars represent the Icr combined statistical and systematic 
uncertainties. The dashed line and stars represent Ling's model and QARM, respectively. References to previous observations are provided 
in the text. 



cause the contribution of cosmic gamma rays that are 
scattered in the atmosphere is not neghgible in the sub- 
MeV/MeV region, we also use Geant4 to simulate the 
transport in the atmosphere as a function of energy, 
zenith angle, and atmospheric depth. Using these simula- 
tions, we analyzed the growth curve in the 125-1250 keV 
range and obtained the fluxes of diffuse cosmic gamma 
rays and atmospheric gamma rays. Our results are con- 
sistent with those of the previous experiments. The flight 
of SMILE-I demonstrated that our ETCC can perform 
gamma-ray selection and powerful background rejection, 



and thus, may improve MeV gamma-ray astronomy by 
having better continuum sensitivity. As the next step, 
we plan to observe the Crab Nebula by using an ETCC 
with a 30 X 30 X 30 cm"^ TPC to demonstrate gamma- 
ray imaging. Subsequently, we plan to observe several 
celestial sources with a long-duration balloon flight, and 
finally aim for an all-sky survey at energies between 150 
keV and 20 MeV with sensitivity up to an order of magni- 
tude better than the current and previous missions, using 
an ETCC with ARM resolution of 2 degrees (FWHM) in 
the MeV region and a detection volume of 50 x 50 x 50 
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Fig. 24. — Continuum detection sensitivity as a function of en- 
ergy with an observation time of 10® s at a significance of 3a. For 
comparison with previous observations, the sensitivities of SMILE 
were overlaid on Fig. 1 of Schonfelder (2004). The sohd and short- 
dashed lines represent the sensitivity of an ETCC with a volume 
of 50 X 50 X 50 cm'^ aboard a long-duration balloon (LDB) through 
Compton scattering, without and with an absorber of Compton- 
recoil electrons, respectively. 



cm"^ onboard a long-duration balloon or satellite. 
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